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QUIESCENT COMPACT GALAXIES AT INTERMEDIATE REDSHIFT IN THE COSMOS FIELD. THE 

NUMBER DENSITY 

IVANA DAMJANOv\ MARGARET J. GeLLER^. H. JABRAN ZAHID^. HO SeONG HwANG^'^ 

ABSTRACT 

We investigate the evolution of compact galaxy number density over the redshift range 0.2 < z < 0.8. 

Our sample consists of galaxies with secure spectroscopic redshifts observed in the COSMOS field. 
With the large uncertainties, the compact galaxy number density trend with redshift is consistent 
with a constant value over the interval 0.2 < z < 0.8. Our number density estimates are similar to the 
estimates at z > 1 for equivalently selected compact samples. Small variations in the abundance of 
the COSMOS compact sources as a function of redshift correspond to known structures in the field. 

The constancy of the compact galaxy number density is robust and insensitive to the compactness 
threshold or the stellar mass range (for M* > 10^° Mq). To maintain constant number density any 
size growth of high-redshift compact systems with decreasing redshift must be balanced by formation 
of quiescent compact systems at z < 1. 

Subject headings: galaxies: evolution — galaxies: fundamental parameters — galaxies: stellar content 
— galaxies: structure 


1. INTRODUCTION 

Structural evolution of quiescent galaxies provides im¬ 
portant insights into the mass assembly of the most mas¬ 
sive galaxies in the local universe. A suite of studies 
have demonstrated that these massive quiescent galaxies 
at z > I are several times smaller than typical similarly 
massive passive systems at z ~ 0 (e.g., Daddi et al. 2005; 
Longhetti et al. 2007; Trujillo et al. 2007; Toft et al. 2007; 
Zirm et al. 2007; Cimatti et al. 2008; van Dokkum et al. 
2008; Buitrago et al. 2008; Damjanov et al. 2009, 2011; 
McLure et al. 2012; van de Sande et al. 2013; Belli et al. 
2014a). Potential sources of systematic errors in both 
galaxy size measurements (based on the surface bright¬ 
ness profile htting) and galaxy stellar mass estimates do 
not signihcantly affect these findings (Muzzin et al. 2009; 
Szomoru et al. 2012). 

Evolution in the abundance of compact systems with 
cosmic time remains an open issue, particularly at 0 < 
z < 1. At high redshift (z > 1), the number density 
of the compact quiescent sources is in the range 10“^ — 
lO^'^Mpc”^ (e.g., Bezanson et al. 2009; Saracco et al. 
2010; Cassata et al. 2011, 2013; Stefanon et al. 2013; 
Barro et al. 2013; van der Wei et al. 2014). The exact 
value of the number density depends on the definition of 
compactness and on the minimal stellar mass probed in 
the study. 

In the local universe, results vary dramatically among 
surveys. Based on the Sloan Digital Sky Survey (SDSS), 
the number of massive compact systems at z ~ 0 is neg¬ 
ligible; the number density is more than three orders of 
magnitude below the high redshift values (Trujillo et al. 
2009; Taylor et al. 2010). The apparent disappearance 
of the massive compact systems from z > 1 to z = 0 
requires a set of dry mergers or alternative exotic mech- 
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anisms to drive galaxy size growth (e.g., Nipoti et al. 
2009a,b; Oogi & Habe 2013; Belli et al. 2014b). How¬ 
ever, the apparent sharp decline in compact galaxy num¬ 
ber density may be the result of biased galaxy selection 
and poor imaging quality of the SDSS (Taylor et al. 2010; 
Carollo et al. 2013). 

In contrast with the SDSS-based studies, results from 
the PM2GC survey at z < 0.1 (Calvi et al. 2011) suggest 
that the abundance of local compact systems is compara¬ 
ble with the number density at high redshift (Poggianti 
et al. 2013b, a). The WINGS survey of nearby clusters 
(Fasano et al. 2006) shows high fractions of compact sys¬ 
tems in galaxy cluster environments (Valentinuzzi et al. 
2010), suggesting a link between compact massive sys¬ 
tems and the galaxy local density. 

Using a sample of point-source objects from the SDSS 
photometric database with signatures of passive evolu¬ 
tion in their Baryon Oscillation Spectroscopic Survey 
spectra (BOSS, Dawson et al. 2013) we identified a large 
sample of intermediate-redshift massive compact galax¬ 
ies (Damjanov et al. 2013, 2014). These objects are at 
least two times smaller than their massive z ^ 0 SDSS 
analogs. Hard lower limits on the number densities of 
these systems at 0.2 < z < 0.6 clearly show that the 
abundance of massive compact systems does not decline 
steeply in this redshift range (Damjanov et al. 2014). 

The intermediate redshift regime is an important link 
between compact samples in the nearby and high-redshift 
universe. The range 0.2 < z < 0.8 presents a sweet spot 
for identifying the compact quiescent galaxy population 
and for probing their structural, dynamical and environ¬ 
mental properties. The volumes surveyed in this red¬ 
shift range are large and compact massive systems are 
bright enough to be targeted by large spectroscopic sur¬ 
veys. Spectro-photometric studies of the intermediate- 
redshift compact population can further constrain the 
mechanisms governing the compact massive galaxy for¬ 
mation and evolution. 

Here we examine a sample of massive (M* > 10^° Mq) 
quiescent galaxies in the COSMOS field with available 
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size measurements based on surface brightness profile fit¬ 
ting and with known spectroscopic redshifts to measure 
the compact galaxy number density at 0.2 < z < 0.8. 
In Section 2 we introduce the parent sample of quiescent 
galaxies and describe the selection functions we use to ex¬ 
tract compact galaxy samples. In Section 3 we select one 
compact sample to demonstrate the method we employ 
to correct for observational selection effects and to esti¬ 
mate the compact galaxy number density. In Section 4 
we examine the effects that the cosmic variance and 
compact galaxy definition have on derived intermediate- 
redshift number densities. In an accompanying paper 
we investigate the Fundamental Plane defined by the 
most massive intermediate-redshift compact systems in 
the COSMOS field (Zahid et al. 2015, , hereafter Pa¬ 
per II). We adopt a cosmological model with Oa = 0.7, 
flM = 0.3, and Hq = 70 km s“^ Mpc“^. 

2. IDENTIFYING MASSIVE COMPACT GALAXIES IN THE 
COSMOS FIELD 

The public data for the COSMOS field offer a unique 
opportunity to identify massive compact galaxies and 
to explore the evolution in their number density as a 
function of redshift. The available datasets include: 1) 
the catalog of structural parameters for objects in the 
1.64 deg^ field covered by the HST ACS F814W imag¬ 
ing (Scarlata et al. 2007; Sargent et al. 2007; Koekemoer 
et al. 2007)^, 2) photometric catalog with stellar masses 
for galaxies in the 1.5 deg^ (unmasked) held of the Ul- 
traVISTA survey (Ilbert et al. 2013; McCracken et al. 
2012)®, and 3) a compilation of available spectroscopic 
redshifts in the COSMOS held (Davies et al. 2014)®. 

Here we describe our approach to constructing a cat¬ 
alog of compact massive galaxies in the redshift range 
0.2< 2 < 0.8. Because the identihcation may depend on 
the approach to size measurement, we remeasure sizes 
for a subset of the smallest objects. We then extract 
catalogs based on the various dehnitions of compactness 
and we summarize the samples. 

2.1. Sample Selection 

We select the parent quiescent sample using the fol¬ 
lowing criteria: 

• Galaxies are in the redshift range 0.2 ^ z ^ 0.8; 

• Galaxy stellar masses are M* > 10^® Mq; 

• Galaxies are dehned as passive: {NUV — r) > 3 x 
(r — J)-1-1 and {NUV — r) > 3.1 (Ilbert et al. 2010, 
2013). 

In addition, we select only intermediate redshift quies¬ 
cent galaxies with robust spectroscopic redshifts reported 
in the Davies et al. (2014) redshift catalog. These objects 
include: 

• zCOSMOS-observed sources with redshifts classi- 
hed as 99% reliable (based on the classihcation 
scheme described in Lilly et al. 2009), 

http: //irsa. ipac. caltech. edu/data/COSMOS/tables/ 
morphology/cosmos_morph_zurlch_l. 0. tbl 
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• Targets with AUTOZ redshifts that are confirmed 
by PRIsm MUlti-Object Survey (PRIMUS, Coil 
et al. 2011; Cool et al. 2013) or zCOSMOS sur¬ 
veys (parameter Z_GEN= 1, 2 in the Davies et al. 
2014 catalog), 

• VIMOS VLT Deep Survey (Le Fevre et al. 2005) 
targets, 

• SDSS DRIO (Ahn et al. 2014) spectroscopic 
sources. 

We crossmatch this sample with the COSMOS catalog 
of galaxy morphology based on Gim2D (Simard et al. 
2002) fitting, using the cutoff magnitude /uvista(AR) = 
23 mag to include only systems with reliable size mea¬ 
surements (Sargent et al. 2007). The final parent sample 
contains 1599 massive quiescent galaxies with available 
size measurements and with spectroscopic redshifts in 
the range 0.2 ^ z ^ 0.8. 

Figure 1 shows the distribution of the COSMOS par¬ 
ent galaxy sample properties. The spectroscopic red¬ 
shift distribution in Figure la shows prominent peaks in 
galaxy number counts at z ^ 0.35, 0.65 < z < 0.7, and 
z ~ 0.725, where there are significant galaxy overdensi¬ 
ties (e.g., Scoville et al. 2007; Guzzo et al. 2007; Kovac 
et al. 2010; Masters et al. 2011; Scoville et al. 2013). 

2.2. Stellar Mass Estimates 

We adopt the stellar masses given in the Ilbert et al. 
(2013) cataolog. Ilbert et al. (2013) estimate the stel¬ 
lar masses for COSMOS galaxies by comparing stellar 
population synthesis models to 30 bands of ultraviolet 
to infrared photometry. The multi-band SEDs are fit 
using the LePHARE^ code written by Arnouts & Ilbert 
(Arnouts et al. 1999; Ilbert et al. 2006). A library of 
synthetic spectra are generated based on the Bruzual 
& Chariot (2003) stellar population models. Synthetic 
SEDs are reddened using the Calzetti et al. (2000) extinc¬ 
tion law with a range of E{B — V) values. Exponential 
star formation histories are adopted with r = 0 — 30 Gyr 
and models span three metallicities ([0.2, 0.4,1] x Zq). 
The stellar mass is the scale factor that minimizes the dif¬ 
ference between the synthetic and observed SEDs. The 
median of the stellar mass distribution is the stellar mass 
given in the Ilbert et al. (2013) catalog. Eor consis¬ 
tency, we use the LePHARE code with similar param¬ 
eters to calculate the stellar masses of SDSS galaxies in 
Section 3.1. 

Eigure lb shows the histogram of the stellar mass dis¬ 
tribution for the parent sample. The blue and red his¬ 
tograms show the mass distribution for galaxies with red¬ 
shifts below and above the median redshift (z ~ 0.5), 
respectively. The galaxies differ in the expected sense 
that low mass galaxies are relatively less abundant in 
the higher redshift sample. In general the decline in the 
mass distribution below the characteristic mass (where 
the mass function of quiescent population itself stays rel¬ 
atively constant, Ilbert et al. 2013) reflects observational 
selection that we must account for in calculating the to¬ 
tal number density of objects as a function of redshift 
(see Section 3.2). 

^ http://www.cfht.hawaii.edu/-arnouts/LEPHARE/lephare. 
html 
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Figure 1. The distribution of spectroscopic redshifts (a), stellar masses (b) and circularized effective radii (c) for the intermediate redshift 
massive quiescent galaxies selected in COSMOS field. Dashed lines indicate median values for each quantity. The stellar mass and effective 
radius panels show low redshift (0.2 z < 0.5; blue) and high redshift (0.5 ^5 2 : ^ 0.8; red) subsamples. 


2.3. Size Measurements 

Robustness of the size measurements is a crucial el¬ 
ement in understanding the nature of compact massive 
galaxies and their connection to the entire quiescent pop¬ 
ulation. Next we check that the sizes of some of the most 
compact objects in the COSMOS field are insensitive to 
the approach to size determination. 

To test size estimates of Sargent et al. (2007), we 
choose 41 quiescent COSMOS galaxies with measured 
sizes and stellar masses, but which are classified as point 
sources in the SDSS DRIO photometric catalog. We in¬ 
clude only the brightest objects with rsDSS < 21.8 to 
match the selection criteria we previously used to esti¬ 
mate the lower limit on the number density of intermedi¬ 
ate redshift compact systems in BOSS (Damjanov et al. 
2014). 

Using Galfit software v3.0.4 (Peng et al. 2010), we fit 
the HST ACS F814W images of SDSS/COSMOS targets 
with two sets of two-dimensional models: a single Sersic 
profile and a composite model that may include up to 
three Sersic profiles. The functional form of the single¬ 
profile surface brightness model is: 


where is the effective radius encompassing half of the 
total flux, Ee is the surface brightness at Re, n is Sersic 
index describing the central concentration, and K{n) is a 
normalization factor. For each galaxy we select several 
(2 — 5) bright stars within a 45" radius and use IRAF 
daophot routines to construct corresponding PSF. Gal¬ 
fit convolves this PSF with a Sersic surface brightness 
profile and uses the Levenberg-Marquardt algorithm to 
match this analytic model to the observed surface bright¬ 
ness prohle. We compare the resulting best-fit circular¬ 
ized effective radius = Re^/q (where q is the axial 


ratio of the best-fit model) to the same parameter of 
the best-fit Gim2D model, The two quantities 

are in excellent agreement (circles in Figure 2a), with 
^Gim2D/^Galfit _ ^_qq q Thus the two methods 
based on the comparison between the Sersic model and 
observed surface brightness profiles, Gim2D and Galfit, 
provide almost identical size estimates. 

We compare the effective radius along the major axis 
of the best-fit single-Sersic profile with the effective ra¬ 
dius within the elliptical aperture (measured using the 
ZEST+ algorithm, Scarlata et al. 2007) for selected ob¬ 
jects in Figure 2b. These two quantities also agree very 
well: ^ o.94 ± 0.17. However, in the 

regime of high surface brightness and compact sizes, Gal¬ 
fit measurements (or similar techniques based on match¬ 
ing with a PSF-convolved model) are more reliable for 
recovering intrinsic sizes of simulated galaxies (Carollo 
et al. 2013). 

Because the HST ACS images reveal a bulge-l-disk 
structure for most of these targets, we also fit the ob¬ 
served surface brightness prohles of compact massive 
galaxies in our COSMOS/SDSS sample with a compos¬ 
ite models. Indeed, only three out of 41 objects are well- 
fitted with a single bulge. After obtaining the best-ht 
multi-S&sic profile for each target, we use the ellipse task 
in the IRAF stsdas.analysis.isophote package to extract 
the intrinsic one-dimensional profile (i.e., before the PSF 
convolution) and to perform a curve of growth analysis. 

The resulting multi-profile effective radius is on average 
very close to the single-profile model size: the mean value 
of the ratio ijmulti-Gafflt/^single-Galfit ^ q 95 (gq^areS 

in Figure 2a). However, lO(eight) out of 41 objects ex¬ 
hibit multi-profile sizes that are more than 10 %{ 20 %) 
smaller than their single-Sersic model sizes. The compos¬ 
ite structure of these galaxies includes a faint disk compo¬ 
nent surrounding a prominent highly concentrated bulge 
which drives the intrinsic effective radius of the multi- 


E(r) = Eg X exp 
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Figure 2. (a): Comparison between circularized effective radii of the best fit Galfit and Gim2D models for 41 intermediate redshift 

compact galaxies in the COSMOS, (b): Effective radius along the major axis of the best-fit Galfit model versus the ZEST+ semi-major 
axis length of the ellipse encompassing 50% of total light. In both panels circles correspond to two single Sersic (1968) profile models. The 
j/—coordinate for squares is the effective radius of the best-fit multi-Sersic profile based on Galfit. All points are color-coded by redshift. 


Sersic profile down to 37% of the single-Sersic profile size 
(for the most extreme case). In contrast, only 5(1) ob¬ 
jects have multi-profile models which are 10%(20%) more 
extended than their single-Sersic fits. The ratio between 
the multi-profile Galfit effective radius and the ZEST+ 
size shows a similar trend (squares in Figure 2b). These 
observed trends confirm the results based on size mea¬ 
surements of simulated multi-component galaxies: mod¬ 
eling such systems with a single Sersic component pro¬ 
duces slightly overestimated sizes (Davari et al. 2014). 
Although almost all of the sample galaxies show a disk 
component in addition to a prominent bulge, their com¬ 
posite structure is at least as compact as their best-fit 
single-profile model. 

Figure Ic shows the size distribution of the massive qui¬ 
escent galaxies. They span the range 0.3 kpc < ^ 

27 kpc with a median value of i?e,c = 2 kpc. The size dis¬ 
tributions above and below the median redshift are more 
similar than the mass distributions. This similarity sug¬ 
gests that correction for an observational bias related to 
size measurements in this redshift range is unnecessary. 

2.4. The Compact Galaxy Catalogs 

We construct several samples of 0.2 < z < 0.8 compact 
galaxies using different compactness criteria defined in 
the literature. All of these definitions are based on the 
position of selected objects in the effective radius vs. stel¬ 
lar mass diagram. They account for the positive trend 
that galaxy sizes exhibit as a function of stellar mass. 
Table 1 summarizes the properties of each selected sam¬ 
ple. 

The fraction of compact systems in each parent sample 
depends on the adopted definition. In Damjanov et al. 
(2014) we pre-select compact galaxy candidates by look¬ 
ing for r{AB) < 21.8 mag systems that are classified 


as point sources in the SDSS photometric database, but 
that show spectroscopic signatures of redshifted quies¬ 
cent galaxies. All of the candidates with available high- 
resolution imaging are at least two times smaller than 
a typical quiescent SDSS galaxy of equivalent mass at 
z ^ 0. If we use the same criteria to select compact 
galaxies in the quiescent sample of the COSMOS field 
(Section 2.3), almost all of candidates are indeed more 
than two times smaller than their massive SDSS coun¬ 
terparts at z ~ 0 (Definition 1 in Table 1). Interestingly, 
this sample contains mostly lower-mass systems, with an 
average stellar mass of (M*) ~ 2.6 x 10^° M©. Because 
of the stellar mass - size relation, massive galaxies can 
be compact even if they appear extended in poor see¬ 
ing conditions (although their size measurement based 
on such imaging may not be accurate). Thus we extract 
different samples using all of the intermediate-redshift 
quiescent galaxies in the COSMOS field with available 
spectroscopic redshifts, stellar masses, and structural pa¬ 
rameters. 

The upper limit on the mass-normalized galaxy size is 
crucially important in determining the number of galax¬ 
ies in each sample. Cutoffs that are close to the size-mass 
relation defined by the local sample produce subsets that 
contain half of the total quiescent galaxy sample (Defini¬ 
tions 2 and 3 in Table 1, Barro et al. 2013 and compact 
galaxies in Cassata et al. 2011). In contrast, defining 
ultra-compact galaxies according to Cassata et al. (2011) 
(Definition 4 in Table 1) puts these systems 0.4 dex (or 
lower) below similarly massive local systems. Although 
this approach naturally produces a much smaller com¬ 
pact sample, mass distributions of systems classified ac¬ 
cording to Definitions 2, 3, and 4 are remarkably similar 
because all three cutoffs are based on the linear log(size)- 
log(stellar mass) relations with similar slopes. 
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Table 1 

The properties of compact COSMOS samples 


Definition 

Reference 

Sample size^ 

(mq)) 

(4) 

^log Af # 

(Pe)^ 


^min 

^max 

( 1 ) 

( 2 ) 

(3) 

(5) 

[kpc] 

( 6 ) 

[kpc] 

(7) 

( 8 ) 

(9) 

(1) log (i?e.c[kpc]) < 0.568 X log (AL*[Mq])‘^ 
-5.74 

Damjanov et al. (2014) 
This work 

37 {90%)=^ 

10.42 

0.27 

1.13 

0.65 

0.24 

0 . 66 ® 

(2) Si ,5 _ log (^^ 1,5 5 ] j > 10.3 

Barro et al. (2013) 
Poggianti et al. (2013a) 

733 (46%)^ 

10.87 

0.38 

2.07 

1.33 

0.20 

0.80 

(3) log (Re,c [kpc]) < 0.54 X log (Af* [MqJ) 

-5.5 

Cassata et al. (2011) 
Cassata et al. (2013) 

840 (52%)s 

10.83 

0.37 

1.83 

1.02 

0.20 

0.80 

(4) log [Re,c [kpc]) < 0.54 X log (M* [MqJ) 

-5.8 

Cassata et al. (2011) 
Cassata et al. (2013) 

82 (5%)f5 

10.82 

0.38 

1.04 

0.58 

0.22 

0.80 


van der Wei et al. (2014) 

215 (13%)f 

11.09 

0.34 

2.79 

1.61 

0.22 

0.80 

f7i,T Finn ^ 1-5 

van der Wei et al. (2014) 

19 (l%)f 

11.38 

0.27 

2.72 

1.08 

0.22 

0.76 


^ Size of the sample is given both as the total number of compact systems and as the fraction of a parent sample. 

^ Effective radius is circularized in all definitions but the ones used in van der Wei et al. 2014 (Definitions 5 and 6). 

^ This definition is based on the best fit linear relation between the size and the mass of 2 ~ 0 quiescent galaxies in our SDSS comparison sample. Compact 
galaxies are defined to be at least a factor of two smaller than SDSS galaxies of the same mass. 

The total number of galaxies in the parent sample - intermediate-redshift quiescent COSMOS galaxies that are classified as point sources in the SDSS 
photometric database - is 41 (Section 2.3). 

® The cutoff magnitude for this sample of r{AB) — 21.8 mag strongly affects the upper redshift limit of the selected sample. 

^ The total number of galaxies in the parent sample - intermediate-redshift quiescent COSMOS galaxies with spectroscopically confirmed redshifts and 
I{AB) < 23 mag - is 1599 (Section 2.1). 

® The total number of galaxies in the parent sample - intermediate-redshift quiescent COSMOS galaxies with spectroscopically confirmed redshifts and 
I{AB) < 23 mag scaled to Salpeter IMF and updated version of Bruzual & Chariot (2003) models - is 1613. 


The steepest d(logi?e)/(i(logM*) gradient (Definitions 

5 and 6 in Table 1, van der Wei et al. 2014) produces 
a compact sample shifted towards higher stellar masses 
and larger sizes. In selecting compact galaxy samples 
van der Wei et al. (2014) use galaxy size measured along 
the galaxy major axis; they do not circularize it. This 
approach is the major factor accounting for the discrep¬ 
ancy between the linear coefficient of Definitions 5 and 

6 and others listed in Table 1. Although mass distribu¬ 
tions (and subsequently size distributions) vary among 
compact samples listed in Table 1, all samples with the 
same limiting magnitude {I{AB)iim = 23 mag) cover the 
full range of redshifts that we explore. None of the tested 
definitions produces a sample of compact systems which 
disappear at low redshift. 

3. THE ABUNDANCE OF COMPACT MASSIVE GALAXIES 
IN COSMOS 

In order to robustly estimate the number density of 
compact galaxies at intermediate redshift we need: a) 
a large enough survey area to provide numerous objects 
and to minimize the effect of cosmic variance, b) high- 
quality imaging to enable reliable size measurements of 
these extreme systems, and c) spectroscopy for secure 
redshift estimates. The COSMOS field, with HST ACS 
imaging available for 1.65 deg^ and its (sparse) spec¬ 
troscopic coverage, is currently one of the best regions 
for this type of analysis. One of compact samples (Sec¬ 
tion 2.4) serves here to illustrate all of the steps in cal¬ 
culating compact galaxy number densities. We demon¬ 
strate that the trend in number density with redshift is 
consistent with an approximately constant abundance in 
the redshift range 0.2 < z < 0 . 8 . 

3.1. Properties of the compact sample 



10.0 10.5 11.0 11.5 12.0 
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Figure 3. Circularized effective radius of single-profile Sersic 
models as a function of stellar mass for the COSMOS compact 
galaxies (circles color-coded by redshift). The red 2D histogram 
represents the distribution of the parent quiescent sample at inter¬ 
mediate redshifts and the red dashed line shows the compactness 
cutoff. The gray 2D histogram shows a set of z ~ 0 SDSS massive 
quiescent galaxies. 


We use Definition 2 from Table 1 to extract a sam¬ 
ple of compact massive quiescent galaxies based on their 
pseudo-stellar mass surface density £ 1 , 5 . This criterion 
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Figure 4. Rest-frame NUV — r vs. rest-frame r — J color for 
the COSMOS sample of compact galaxies (circles color-coded by 
redshift). The gray 2D histogram show the underlying distribution 
of the parent sample. The black dashed lines show the boundary 
between star-forming and quiescent systems (Ilbert et al. 2010). 


was introduced by Barro et al. (2013) who applied it to 
estimate the number density evolution of compact galax¬ 
ies at 0.5 < z < 3 in the CANDELS survey (Grogin et al. 
2011; Koekemoer et al. 2011). Poggianti et al. (2013a) 
employed the same definition to show that the number 
density of compact massive systems in the nearby Uni¬ 
verse (z < 0.1 galaxies drawn from the PM2GC survey) 
is very similar to their abundance at high redshift. 

Figure 3 shows the position of selected compact sys¬ 
tems (circles colored by redshift) in the size-stellar mass 
parameter space compared with the distribution of all 
massive quiescent COSMOS systems at 0.2. < z < 0.8 
(red two-dimensional (2D) histogram). For the compar¬ 
ison z ~ 0 sample (gray scale 2D histogram) we em¬ 
ploy a set of spectroscopically confirmed quiescent galax¬ 
ies from SDSSS DR7 with available size measurements 
(Damjanov et al. 2014; Simard et al. 2011). We derive 
stellar masses for the z ^ 0 quiescent sample using the 
same spectral energy distribution fitting method as Il¬ 
bert et al. (2013). The COSMOS compact sample con¬ 
tains 46% of the parent quiescent sample at intermedi¬ 
ate redshift (Table 1). These dense intermediate-redshift 
systems are smaller than 95% of local massive passive 
galaxies. 

In rest-frame color-color space the compact sample co¬ 
incides, by design, with the locus of quiescent systems 
(Figure 4, Section 2.1). Compact galaxies at lower red- 
shifts show redder r — J and NUV — r rest-frame col¬ 
ors because they are dominated by older stellar popu¬ 
lations. Both galaxy colors become bluer (i.e., galaxies 
are younger) as their redshift increases (see also Paper 
II). Although this trend is clear (Figure 4), there is also 
some evidence of mixing. A small fraction of compact 
galaxies at the lower redshift end - ~ 5% of all compacts 
at z ^ 0.4 - show bluer colors (with NUV — r < 4.5 and 
r — J < I) and appear to be more recently quenched. On 


the other hand, ~ 12 % of compact systems at z ^ 0.6 
have very red rest-frame colors - NUV — r > 5.5 and 
r - J > 1. 

Figure 5 shows the distribution of redshifts, masses and 
sizes for the compact sample. The median redshift of the 
compact galaxies (z<.ompact = 0.65, Figure 5a) exceeds 
that for the parent sample (zparent = 0.53, Figure la). If 
compact systems are preferentially located in overdense 
regions, this difference may be related to the presence 
of the largest dense structure in the COSMOS held at 
z ^ 0.7 We can test this conjecture only after correcting 
both the parent and compact samples for spectroscopic 
incompleteness. 

The size distributions of the compact sample below and 
above its median redshift are very similar. The median 
compact sizes in two redshift bins are almost identical 
{Re,c ~ 1-7 kpc, Figure 5c); they are naturally smaller 
than the sizes for the parent sample (Figure Ic). 

3.2. Observational selection effects 

The sampling by spectroscopic surveys targeting the 
COSMOS held is sparse at intermediate redshift and we 
must account for this observational incompleteness to de¬ 
rive compact galaxy number densities. In addition, a 
magnitude limited sample with available size measure¬ 
ments (Sargent et al. 2007) can have a mass distribution 
that departs from the true mass function of quiescent sys¬ 
tems in the lowest mass bins (M* ^ I0^° Mq) at higher 
redshifts (z > 0.5). Here we describe the two corrections 
required to calculate the intermediate-redshift number 
densities of massive compact systems. 

To estimate the fraction of compact galaxies with avail¬ 
able spectroscopic redshifts in the underlying quiescent 
population, we use a method based on the position of 
our spectroscopic targets in color-color-magnitude space 
(Damjanov et al. 2014). Here we add galaxy mass as the 
fourth dimension of the parameter space we explore. For 
each galaxy in the compact sample the (in)completeness 
correction factor is the ratio between the number of 
spectro-photometric and the number of photometric ex¬ 
tended sources (i.e., galaxies) with similar colors, appar¬ 
ent magnitude and stellar mass as the selected system. 
This approach is appropriate if the spectroscopic sample 
is a compilation of sparse surveys with different selection 
criteria, such as our sample of COSMOS galaxies with 
spectroscopic redshifts (Section 2.1). 

We construct the four-dimensional (4D) parameter 
space using galaxy rest-frame colors NUV — r and r — J, 
that are also used to discriminate between star-forming 
and quiescent galaxies (Section 2.1), apparent /—band 
magnitude, and stellar mass. For each compact galaxy 
we define a cell with A{NUV — r) = A(r— J) = 0.2 mag, 
Al — 0.4 mag, Alog(M*[M 0 ]) = 0.2 sides centered at 
the position of the compact target in this parameter 
space. We count galaxies with measured sizes located 
within a given cell and divide their number by the num¬ 
ber of galaxies with additional spectroscopic information. 
This ratio is then the estimated spectroscopic complete¬ 
ness correction at the position of selected compact target. 

We illustrate our method in the panels of Figure 6 that 
show compact objects (circles) in the projections of the 
{NUV — r rest-frame color - r—J rest-frame color - I ap¬ 
parent magnitude - stellar mass) parameter space. Each 
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Redshift log(M, [Mq]) log(R, Jkpc]) 


Figure 5. Distribution of spectroscopic redshifts (a), stellar masses (b) and circularized effective radii (c) of the intermediate redshift 
massive quiescent galaxies in COSMOS classified as compact based on their pseudo-stellar mass surface density E 1.5 (Definition 2 in 
Table 1). Dashed lines show median values for each quantity. We show stellar mass and effective radius panels subsamples below (blue) 
and above (red) median redshift (z = 0.65). 
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log(M.IMal) log(M.|M3|) log(M.[Ma]) 



Spectroscopic completeness 

Figure 6 . Positions of compact galaxies (selected using Definition 2 from Table 1) in six projections of the parameter space we use to 
calculate spectroscopic completeness factors. These parameters include: rest-frame NUV — r, observed I magnitude, rest-frame r — J, and 
galaxy stellar mass M*. Circles color-coded by the spectroscopic completeness factor represent compact systems. The open rectangle in 
each panel shows the projected size of the cell that we use to calculate the spectroscopic completeness level. See the text for details. 
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Figure 7. Corrected distribution of stellar masses (a) and cir¬ 
cularized effective radii (b) for a sample of intermediate redshift 
massive quiescent galaxies in the COSMOS field at low redshift 
(0.2 ^ 2 : < 0.5; blue) and high redshift (0.5 ^ 2 : ^ 0.8; red). 
Dashed lines indicate median values for each quantity. See the 
text for details. 



Figure 8. Mass function of the corrected parent quiescent sample 
based on the COSMOS spectroscopic dataset (points) compared to 
the mass function defined by the COSMOS quiescent systems with 
photometric redshifts (Ilbert et al. 2013, , shaded regions), (a) 
and (b) panels correspond to the low- 2 : and high- 2 ; subsamples, 
respectively. 

circle is color-coded by the spectroscopic correction we 
apply to the corresponding object. The open rectangles 
depict the size of the 4D cell in all six projections. In 
practice we construct a separate 4D cell around the po¬ 
sition of each compact target defined by its stellar mass, 
rest-frame colors, and /—band magnitude. 

The size of the cells we use to estimate the correc¬ 
tion for spectroscopic completeness is set by the num¬ 
ber and mass distribution of galaxies in the photometric 
sample. Using a grid of cell sizes, we apply a range of 
corrections for spectroscopic incompleteness to the par¬ 
ent sample of quiescent galaxies with measured sizes. As 
a result, the same galaxy has different incompleteness 
factors associated with it in different corrected parent 
samples. We divide each corrected parent sample in two 
subsets defined by spectroscopic redshift (0.2 < z < 0.5 
and 0.5 ^ z ^ 0.8). For both redshift intervals we re¬ 
quire the corrected number of galaxies in each narrow 
mass bin (A log(M* [Mq]) = 0.2) to be close to (but not 
to exceed) the number of similarly massive galaxies with 
photometric redshifts in the same broad redshift range. 
This procedure selects the optimal cell size which we then 
use to derive correction factor for each system in our 
compact COSMOS sample. The resulting spectroscopic 
completeness factors cover the range 0.1 — 1 (as shown in 
Figure 6) with an average value (and standard deviation) 


of 0.57 ±0.26. 

In addition to its spectroscopic incompleteness, our 
compact sample may be incomplete at the lower mass end 
because of the magnitude limit of the catalog containing 
the galaxy structural parameters (lumiAB) = 23). We 
account for this additional incompleteness by dividing 
galaxies from the COSMOS photometric catalog (Ilbert 
et al. 2013) into low-redshift (0.2 < Zphot < 0.5) and 
high-redshift (0.5 < Zphot ^ 0.8) subsets and calculating 
the fraction of galaxies with I{AB) < 23 in the under¬ 
lying quiescent galaxy population in narrow mass bins 
(Alog(M*[MQ]) = 0.2) for both subsets. This ratio is 
less than unity only in two mass bins of the high-z sub¬ 
sample: it is 0.68 at M* 1.2 x 10 ^° Mq and 0.89 at 
M* ~ 2 X 10^° Mq. To test the method we apply the 
correction to the parent galaxy sample, which we first 
correct for the spectroscopic incompleteness. 

Figure 7 shows the mass and size distributions of the 
parent sample of quiescent COSMOS galaxies with spec¬ 
troscopically confirmed redshifts, corrected for the in¬ 
completeness due both to spectroscopic sampling and 
magnitude limit of the morphological catalog. Although 
the shape of the size distribution (Figure 7b) is the same 
as for the input sample (shown in Figure 1), the galaxy 
mass distribution of the corrected sample is flatter in the 
low-mass regime (10^° Mq < M* < 5 x 10^° Mq), at 
least for systems in the low-z subsample (blue histogram 
in Figure 7a). 

The number of galaxies of the high-z galaxy subset 
still declines in the lowest mass bins (red histogram in 
Figure 7a). To check whether the corrected mass func¬ 
tion of our spectroscopic parent sample corresponds to 
the independently observed mass function of COSMOS 
quiescent galaxies spanning the same redshift range (Il¬ 
bert et al. 2013), we compare the shapes of two distri¬ 
butions in Figure 8. In the low-z regime our corrected 
parent spectroscopic sample produces a mass function 
(blue circles) of the same shape as the mass function of 
the COSMOS photo-z sample (blue shaded region). In 
the high-z regime, the mass function of our corrected 
subsample departs somewhat from the shape defined by 
the COSMOS quiescent photo-z sample only in the low¬ 
est mass bin (M*(0.5 ^ z << 0.8) ~ 10^° Mq). At all 
stellar masses M* > 2 x 10^° Mq our spectroscopic and 
magnitude related completeness factors provide a cor¬ 
rected parent sample of 0.2 < z < 0.8 quiescent systems 
with mass functions that agree extremely well with the 
observed mass functions based on an independent, de¬ 
tailed analysis of the COSMOS quiescent galaxies with 
photometric redshifts. 

Finally, we test the corrections applied to the spec¬ 
troscopic parent sample by comparing number densities 
based on corrected galaxy counts in Az = 0.2 redshift 
bins within the 0.2 < z < 0.8 redshift interval with the 
number density evolution of massive quiescent galaxies 
from the COSMOS photometric sample (Carollo et al. 
2013, right panel of their Figure 6). We find excellent 
agreement, with a positive ~ 30% offsets in the redshift 
bin (centered at z = 0.7) covering the largest known over¬ 
density in the COSMOS field, for which Carollo et al. 
(2013) make a correction. This test confirms that the 
spectroscopic and magnitude related completeness fac¬ 
tors do not overestimate the absolute number of COS- 
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MOS galaxies in the corrected spectroscopic parent sam¬ 
ple of massive quiescent systems in the 0.2 < z < 0.8 
redshift range. 

3.3. The number density of compact galaxies 

We apply the correction factors from Section 3.2 to 
the spectroscopic subsample of compact intermediate- 
redshift galaxies in COSMOS selected according to Defi¬ 
nition 2 (Table 1). We then divide the corrected number 
of compact systems into narrow redshift bins (Az = 0.1) 
and compute the corresponding comoving volume ele¬ 
ment to estimate the abundance of compact galaxies as 
a function of redshift in the range 0.2 ^ z ^ 0 . 8 . 

The uncertainty in the number densities of compact 
COSMOS galaxies includes errors in the estimated stel¬ 
lar mass, errors in the size measurement, and the con¬ 
fidence limits from Poisson statistics. We build a set of 
Monte Carlo simulations where we vary both the stel¬ 
lar mass and the size of our selected objects taking 
into account average (Icr) errors of A(log(M*[M 0 ])) = 
A(log(i?e(,c)[kpc])) = O.ldex. The error in stellar mass 
is dominated by the uncertainty due to the choice of stel¬ 
lar library (Santini et al. 2015). The error in size corre¬ 
sponds to the maximum difference between the single- 
Sersic effective radius we measure using Galfit and the 
corresponding size given in the GOSMOS morphological 
catalog (Sargent et al. 2007, see Figure 2). The set of 
1000 realizations provides a range of number densities 
for each redshift bin and lower stellar mass limit. We 
fold standard deviations of the resulting number den¬ 
sity distributions into the total uncertainty estimates to¬ 
gether with the Poisson confidence limits corresponding 
to the Gaussian Icr errors (Gehrels 1986). The hnal er¬ 
ror budget is dominated by the confidence limits based 
on Poisson statistics. 

We compare the resulting evolution of compact galaxy 
number density from z = 0.8 to z = 0.2 with the val¬ 
ues from Poggianti et al. (2013a) at z ~ 0 and Barro 
et al. (2013) at z > 1 using the same compactness 
criterion (Figure 9). At the same stellar mass limit 
(M* > 10^° Mq, blue symbols in Figure 9) and with 
the large uncertainty, the abundance of compact sys¬ 
tems increases by 0.2 dex from z ^ 0 to z ~ 1.5. 
The observed number densities are also consistent with 
a roughly constant value of ~ x 10 “^ Mpc“^ in 

this redshift range (see Section 4.2 and Table 2). Thus 
the number density of compacts in the local universe is 
very similar to the value reported at high redshift (Pog¬ 
gianti et al. 2013a). The largest deviations from this 
constant trend (and the largest absolute values of com¬ 
pact galaxy number densities) occur in the redshift bins 
with the largest of overdensities in the GOSMOS field (at 
z ~ 0.35,0.65,0.75, see Section 2.1). 

To further probe the number density evolution of the 
low-mass and high-mass compact systems, we divide the 
sample into two subsets based on the characteristic stel¬ 
lar mass of the COSMOS mass function for intermediate- 
redshift quiescent galaxies {M* = 8 x 10^° Mq, llbert 
et al. 2013). Figure 10 shows that the number densities of 
low-mass and high-mass compact galaxies exhibit almost 
identical patterns, with absolute number densities of low- 
mass compacts being slightly exceeding (but within Icr 
uncertainty intervals of) the abundances of high-mass 


compact systems in all redshift bins. The only redshift 
range where the two number densities differ more than 
Icr is 0.7 < z ^ 0.8. In this redshift bin the abundance of 
compact galaxies with stellar masses above the charac¬ 
teristic mass M* is ~ 2 times lower than the abundance 
of similarly compact but less massive systems. Interest¬ 
ingly, this redshift range includes the largest structure 
known in the COSMOS field (Scoville et al. 2007, 2013). 

Overall, the number density of compact galaxies se¬ 
lected based on their pseudo-stellar mass surface den¬ 
sity El 5 (Definition 2 in Table 1) is approximately con¬ 
stant for 0.2 < z < 0.8. The study of the ZEST+- 
measured size distribution of COSMOS galaxies in this 
redshift range yields similar results: in the stellar mass 
range covered by the bulk of quiescent galaxy population 
(10.5 < log(M*/M 0 ) < 11 ) the abundance of the dens¬ 
est systems does not evolve (Carollo et al. 2013). The 
abundances of compact systems in the COSMOS field 
at these redshifts are similar to the values reported for 
equivalently defined compact galaxies at both lower and 
higher redshifts. The variations in the number density 
values within this redshift range correspond to the red¬ 
shifts of known COSMOS overdensities. The abundance 
of compact galaxies is approximately constant over this 
redshift range regardless of the galaxy mass in the range 
we explore. 

4. DISCUSSION 

There are several observational effects that can alter 
the size and properties of the selected compact galaxy 
sample. These issues include the underlying survey size 
and the selection criteria used to define the sample. For 
a survey like COSMOS, cosmic variance may produce 
large deviations from the typical value in the number of 
(compact) galaxies in any redshift range. The definition 
of compactness influences the total number of galaxies 
in the sample and their mass distribution. Here we in¬ 
vestigate these two observational issues and their effect 
on the resulting number densities of compact systems in 
COSMOS. 

4.1. Cosmie varianee 

Enhanced number densities of compact galaxies in 
some intermediate-redshift bins correspond to dense 
structures in COSMOS (Section 3.3). Here we compare 
lower limits to the number density of compacts deter¬ 
mined from a very large volume survey, BOSS, to the 
similarly defined limits from COSMOS in order to ex¬ 
amine the potential impact of cosmic variance on the 
resulting number density of compacts. 

Using the approach from Damjanov et al. (2014) we 
select COSMOS sources with r{AB) < 21.8 that ap¬ 
pear as stars in the SDSS images and have spectral en¬ 
ergy distributions (SEDs) of intermediate-redshift mas¬ 
sive quiescent galaxies. In addition, these objects are ex¬ 
tended in the HST ACS images of the COSMOS field 
and have measured sizes. The final SDSS/COSMOS 
compact sample contains a subset of these objects that 
are at least two times smaller than a typical local qui¬ 
escent galaxy of equivalent stellar mass (Definition 1 
from Table 1). In comparison with other samples listed 
in Table 1, the SDSS/COSMOS compact sample has a 
relatively small number of systems with stellar masses 
10^° M 0 < M* < 10^^ M 0 . We use this compact 
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Figure 9. Number density of M, > Mq compact galaxies with Si s > 10.3 Mq kpc“^ ® (Definition 2 in Table 1) as a function of 
redshift for the COSMOS intermediate-redshift compacts (circles), a set of galaxies at 2 : < 0.1 (Poggianti et al. 2013a, diamonds, PM2GC 
sample, ), and for high-redshift CANDELS targets (squares, Barro et al. 2013). The blue line and the blue shaded area correspond to the 
average number density (itlcr) of compact galaxies at 0 < 2 < 2. 


dataset to derive number densities that we label as lower 
limits in order to distinguish SDSS/COSMOS selection 
from the complete sample of compact massive galax¬ 
ies in COSMOS. The number densities of the compact 
SDSS/COSMOS systems are comparable with measured 
number densities of intermediate-redshift compact sys¬ 
tems based on SDSS/BOSS sample, that are also inter¬ 
preted as lower limits because of the selection effects as¬ 
sociated with the parent BOSS sample (Damjanov et al. 
2014). 

Figure 11 displays the excellent agreement between the 
two sets of lower limits on compact galaxy number densi¬ 
ties. The abundances derived from the SDSS/COSMOS 
sample (circles) are above and mostly within the uncer¬ 
tainties of the SDSS/BOSS values (stars) for all stellar 
mass limits (color-coding) in the 0.2 < z < 0.5 redshift 
range. Because the BOSS survey covers an area almost 


4000 times larger than COSMOS and gives similar re¬ 
sults for a similarly selected sample, we conclude that 
the effects of cosmic variance are not dominant. Fur¬ 
thermore, this result shows that the selection algorithm 
of the BOSS quasar survey (Ross et al. 2012), where we 
find compact SDSS/BOSS systems, does not exclude any 
subset of the compact population at these redshifts. 

Interestingly, at some stellar mass limits in this redshift 
range we do not detect any COSMOS compact galaxies 
classified as stars in the SDSS imaging (upper limits cor¬ 
responding to the 2(7 confidence level for non-detection in 
Figure 11). Intermediate-redshift massive galaxies that 
are compact enough to be misclassified as point sources 
in low-resolution images (FWHM(SDSS)^ 1"5) are very 
rare. The search for such extreme objects requires a large 
volume, like the BOSS survey. 

At the high end of the redshift range covered by the 
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Figure 10. Number density of compact galaxies with E 1.5 > 
10.3 Mq kpc“^'^ (Definition 2 in Table 1) as a function of redshift 
for the COSMOS sample. The two subsamples are based on galaxy 
stellar mass. The upper mass limit of the low-mass subset (blue 
circles) and the lower mass limit of the high-mass subset (red cir¬ 
cles) correspond to the characteristic mass of the COSMOS mass 
function for quiescent galaxies at 0.2 < 2 < 0.8 (Ilbert et al. 2013). 

SDSS/COSMOS sample (0.6 < z < 0.8) galaxies become 
too faint for our SDSS apparent magnitude cutoff. Upper 
values on the number densities of SDSS/COSMOS com¬ 
pact systems in this redshift regime reflect the limitation 
of the bright magnitude-limited sample. 

The agreement between the lower limits on the num¬ 
ber densities based on SDSS/COSMOS and SDSS/BOSS 
samples shows that cosmic variance does not dominate 
the resulting overall trend in the COSMOS compact 
galaxy number density with redshift. It also demon¬ 
strates that the broad color selection of the BOSS quasar 
survey is not biased against compact massive intermedi¬ 
ate redshift galaxies. This comparison is a reminder that 
in the intermediate-redshift regime, massive systems may 
appear extended even in low-quality imaging surveys and 
still be defined as compact. In order to select all compact 
quiescent galaxies using a criterion based on the size- 
stellar mass relation, the parent sample has to include 
all massive systems with available size measurements. 

4.2. Compact galaxy definition 

Compact galaxy samples are defined based on the 
galaxy position in the size-stellar mass parameter space 
(Section 2.4). The two major factors affecting the size 
and properties of compact systems are: 1) the slope of 
the best-fit linear relation between galaxy size and stellar 
mass at z ~ 0 (linear coefhcient of the selection function), 
and 2) the distance between the threshold for compact¬ 
ness and the z ^ 0 size-mass relation (zero-point of the 
selection function). Although the slope of the local size- 
mass relation for quiescent systems is well known (~ 0.6; 
e.g., Shen et al. 2003), even a small change in value may 


Figure 11. The number density of compact galaxies as a func¬ 
tion of redshift for two bright {r{AB) < 21.8) intermediate-redshift 
samples: 1) SDSS/BOSS galaxies that appear as point sources 
in the photometric SDSS database and show signatures of red- 
shifted evolved systems in their BOSS spectra (Damjanov et al. 
2014, stars) and 2) COSMOS quiescent galaxies classified as point 
sources in the photometric SDSS database (circles). Each sample is 
represented by a set of five subsets color-coded based on the lower 
mass limit used to select the subset. All upper limits correspond 
to the Poisson upper limit for the 2cr confidence level in the case 
of non-detection. 

produce a very different mass distribution for the com¬ 
pact quiescent sample. The position of the compactness 
threshold crucially affects the number of galaxies in the 
selected sample. We investigate the effects of these two 
factors on the redshift evolution of compact galaxy num¬ 
ber densities for the COSMOS samples (Table 1). 

The absolute value of the compact galaxy abundance 
depends on the number of objects in the sample (i.e, on 
the zero point of the selection function). We perform 
the procedure described in Section 3 using four different 
criteria for identifying compact galaxies in the COSMOS 
sample (Definitions 3-6 in Table 1). Compact COSMOS 
galaxy subsets corresponding to a more extreme cutoff 
for compactness (circles in panels (4) and (6) of Fig¬ 
ure 12) produce an order of magnitude lower number 
densities than samples based on more inclusive selection 
criteria for the same size-mass relation slope (circles in 
panels (3) and (5) of Figure 12, respectively). 

The linear coefhcient of the selection function pro¬ 
duces the spread between number densities of subsets 
with different lower mass limits (color-coded circles in 
Hgure 12). A shallower slope of the size-mass relation 
gives a range of abundances that are progressively lower 
with increasing minimum stellar mass of galaxies in the 
sample. The number density of compact systems with 
M* > 10^° Mq is a factor of four (0.6 dex) larger than 
the number density of the most massive compact galaxies 
with M* > 10^^ Mq (color-coded circles in Figure 9 and 
in panels (3) and (4) of Figure 12). A selection function 
with a lower linear coefficient is more likely to include 
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Figure 12. Number density of compact galaxies as a function of redshift for different compactness thresholds (numbered as in Table 1). 
Each panel contains compact samples defined using the expression shown in that panel. Our results (circles) are clearly similar to the 
results of other surveys that define compact massive galaxies in the same manner (triangles in the upper panels (Cassata et al. 2013) and 
diamonds in the lower panels (van der Wei et al. 2014)), with the exception of the lowest redshift point of van der Wei et al. (2014) in 
panel (5). All results are color-coded by the lower limit on the stellar mass of galaxies in the sample. 


Table 2 

The evolution of compact galaxy number density 


Definition^ 

M, 

0.2 < ^ < 0.8*’’'^ 

0 < 2 < 1.5'"’'* 


O 

o 

1 



(1) 

(2) 

(3) 

(4) 


(2)Si.5_w(4*a [kp"?a]) >10.3 

> 1 

n(z) 

(3) log (i?e,c [kpc]) < 0.54 X log (M, [Mq]) 

> 1 

n(z) 

(4) log (i?e,c [kpc]) < 0.54 X log (M, [Mq]) 

(5) , r ,X0 76 < 2.5 

(A< (lOllM0])“."° 

> 1 

n{z) 

> 5 

n(z) 

(6) , r ,X0 76 < 1.3 

(A< (iO11M0])“."° 

> 5 

n{z) 


10 {- 3 . 55 ± 0 . 26 ) (1 ^ 2 )( 1 - 25 ± 1 - 23 ) n(z} 

10 {- 3 . 34 ± 0 . 26 ) (1 + 2 )( 0 -«± 1 - 25 ) ri(z) 

10 {- 5 . 61 ± 0 . 27 ) (1 + 2 )( 6 . 36 ± 1 . 22 ) 
10 {- 4 . 62 ± 0 . 29 ) (1 + 2 )( 3 - 16 ± 1 . 55 ) n{z) 
10 (- 5 . 46 ± 0 . 42 )(i + ^)( 2 . 26 ± 1 . 90 ) 


]^g{-3.43±0.23) ^ ^'^(0.57±1.07) 

]^q{-3.20±0.16) ^ ,2^(-0.27±0.63) 

]^q{-4.72±0.28)(]^ ^ 2^(2.82±1.12) 
]^q{-4.02±0.19)(]^ j_ 2)(0.07±0.67) 

]^q{-5.04±0.30)(']^ 2)(0.04±1.02) 


^ Equivalent to definitions in Table 1 
Based on COSMOS data only. 

^ Based on the COSMOS field number densities combined with equivalently derived number densities from the literature. 

^ The best-fit parameters and their Irr uncertainties are formal results of the weighted least square linear regression: log(n) 


the maximum fraction of objects at the lower end of the 
parent mass distribution. 

A steeper d{\og Re)/d{\og M^) gradient of the selection 
function excludes low-mass systems and produces a com¬ 
pact sample shifted to higher stellar mass (Section 2.4). 
The resulting number densities for different limiting stel¬ 
lar masses thus exhibit much smaller scatter as a function 
of the difference between the complete compact sample 
(M* > 10^° Mq) and the most massive compact systems, 
in the range of 0 — 0.3 dex (panels (5) and (6) of Fig- 


= a + P log(l -h z) 

ure 12). This difference becomes more prominent in the 
high-redshift bins that probe lager comoving volumes; 
these samples are more likely to include rare extremely 
compact low-mass systems (see panel (6) of Figure 12 for 
the example of selection function with high linear coeffi¬ 
cient and low zero-point). 

We compare derived number densities of these differ¬ 
ently selected samples of compact galaxies in COSMOS 
with the results of other studies that employ equivalent 
selection criteria. Definitions 3 and 4 were introduced by 
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Cassata et al. (2011) and we add their updated results for 
GOODS-South field (Gassata et al. 2013) to panels (3) 
and (4) of Figure 12 (blue triangles). For both definitions 
we find very good agreement between the abundances of 
compact systems in these two fields. In panels (5) and 
( 6 ) of Figure 12 we compare GOSMOS results with mea¬ 
surements based on the sample of equivalently defined 
compact galaxies in the GANDELS fields (van der Wei 
et al. 2014, yellow diamonds). Here we find agreement 
for the abundances of the most extreme (and thus rare) 
compacts (panel ( 6 )) that also have the highest associ¬ 
ated uncertainties. For the more inclusive compactness 
threshold (panel (5)) van der Wei et al. (2014) suggest 
that the number density of compact systems in GAN¬ 
DELS declines drastically between z = 1 and z = 0. As 
we discuss below, the range of possible number density 
trends with redshift for the COSMOS 0.2 < z < 0.8 
dataset in panel (5) of Figure 12 is sufficiently large to 
include the lowest rates of evolution in compact galaxy 
number density for five CANDELS fields reported by 
van der Wei et al. (2014). However, in combination with 
the z > 1 CANDELS data points, the COSMOS field 
data suggest more moderate change in the number den¬ 
sity of compacts in the range 0.2 < z < 1.5. 

To test wether the observed number density evolution 
is consistent with constant abundance of compact galax¬ 
ies in the intermediate redshift range, we derive power- 
law relations between number density and redshift for 
every compact sample in Figures 9 and 12. Taking the 
uncertainties in number densities into account, we per¬ 
form weighted least square fits 1) to the COSMOS data 
in the 0.2 ^ z ^ 0.8 range alone, and 2 ) to combined data 
points from COSMOS and other fields covering redshift 
range 0 < z ^ 1.5. The formal best-fit parameters with 
associated la errors are listed in Table 2. Some formal 
± 1(7 intervals for the (1 ± z) exponent include negative 
values, mostly driven by the highest redshift bin we ex¬ 
plore (z ~ 1.5). Several studies have found that compact 
galaxy number density starts declining with increasing 
redshift somewhere between z ~ 1.5 and z ~ 2 (e.g., 
Barro et al. 2013; van der Wei et al. 2014). 

The number density increases strongly in the 0.2 < 
z < 0.8 redshift range only for the sample selected using 
Definition 4 (ultra-compact galaxies from Cassata et al. 
2011, 2013). This sample contains a small number of 
0.2 < z < 0.8 COSMOS objects (82, see Table 1), most 
of which are in the two highest redshift bins (z ^ 0.65 
and z ~ 0.75) with largest volumes and the densest struc¬ 
tures. Thus the trend in number density with redshift is 
driven by the two high-redshift points with smallest as¬ 
sociated uncertainties. A compilation of COSMOS data 
points and the Cassata et al. (2013) results for ultra¬ 
compact galaxies shows a much less pronounced posi¬ 
tive trend in galaxy number density with redshift. For 
all other definitions of galaxy compactness and in both 
0.2 < z < 0.8 and 0 < z < 1.5 redshift intervals, con¬ 
stancy of the compact galaxy number density (/3 ~ 0) is 
included in the 95% confidence interval for the best-fit 
slope of the linear log(n) oc /31og(l ± z) relation. Fits 
to the data demonstrate that our primary conclusion is 
independent of definition and sampling; the number den¬ 
sity of compact galaxies does not signihcantly evolve in 
the 0 < z < 1.5 redshift range. 

The lack of evolution in number density with redshift 


suggests that either 1 ) compact massive galaxies form 
in the first 4 — 5 Gyrs of cosmic history (down to red¬ 
shift z ~ 1.5) and do not evolve in size after that time 
(static equilibrium), or 2 ) the size growth of high-redshift 
compact galaxies is balanced by the emergence of new 
compact massive galaxies at redshift z < 1 (e.g., Carollo 
et al. 2013, Paper H). To explain the evolution in the av¬ 
erage size of massive galaxy with redshift, many studies 
have suggested that massive galaxies grow in the red¬ 
shift range 0 < z < 1, mainly through a series of minor 
mergers (e.g., Newman et al. 2012; Lopez-Sanjuan et al. 
2012). Thus these objects are obviously removed from 
the compact sample and must be replaced by some mech¬ 
anism. Interestingly, the most massive compact galaxies 
in the SDSS at z ~ 0 are recently quenched systems 
(with ages < 4 Gyr, or a formation redshift Zform ^ 0.6, 
Trujillo et al. 2009; Ferre-Mateu et al. 2012). Further¬ 
more, at least a fraction of compact galaxies at inter¬ 
mediate redshift have E + A-type spectra (i.e, show fea¬ 
tures of young passive systems, Damjanov et al. 2014). 
Thus compact massive galaxies are forming at interme¬ 
diate redshift. However, the fraction of newly quenched 
galaxies in the underlying compact population in this 
redshift range remains an open question. Possible pro¬ 
genitors of intermediate-redshift compact E -\- A systems 
may be compact starbursts evolving in analogy to the 
formation and evolution of compact massive quiescent 
galaxies at high redshift (see e.g.. Toft et al. 2014; Barro 
et al. 2014; Williams et al. 2014; Stefanon et al. 2013; 
Nelson et al. 2014 for observational perspective and e.g., 
Wellons et al. 2015; Zolotov et al. 2014 for theoretical 
consideration). 

5. CONCLUSION 

We use the COSMOS data to measure the abundance 
of massive compact quiescent galaxies at intermediate 
redshift. We correct the parent spectroscopic sample of 
quiescent systems for incompleteness due to both sparse 
spectroscopic sampling and magnitude-limited size mea¬ 
surements to recover the known mass function for quies¬ 
cent galaxies in COSMOS. The trend in number density 
of compact galaxies with redshift in the 0.2 ^ z ^ 0.8 
range is, within uncertainties, consistent with constant 
abundance of these systems at intermediate redshift. 

Systematic observational issues have negligible effect 
on the abundance of compact systems: 

• The dehnition of compact galaxies does not signif¬ 
icantly affect the dependence of their number den¬ 
sity on redshift. 

• Within the uncertainties, the number density of 
compact galaxies does not evolve between z ~ 1.5 
and z ~ 0 . 

• Comparison of a compact sample from BOSS with 
a similar sample based on COSMOS demonstrates 
that cosmic variance does not dominate the results. 

The evolution of compact galaxies at intermediate red¬ 
shift produces a population of objects where 

• The number density of objects does not depend 
strongly on the redshift. 
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• The abundance of compact systems with masses 
M* > 10^° Mq at intermediate redshift ranges 
from 7 X 10“® Mpc“^ to 5 x 10“"^ Mpc“^, depending 
on the criteria for compactness. 

• The trend in compact galaxy number density with 
redshift is insensitive to the adopted minimum stel¬ 
lar mass over the range 10^° Mq < M* < 6 x 
IQI^Mq. 

• Thus, if compacts are indeed the seed for large ob¬ 
jects, some compacts must form at redshift z < 1 
to replenish the supply. 

Analysis of the COSMOS field further suggests that 
the abundance of intermediate-redshift compact objects 
may depend on environment (as indicated by Valentin- 
uzzi et al. 2010 for the z ^ 0 regime and by Raichoor 
et al. 2012 at high redshift). In Paper II we show that 
the most massive compact objects in this redshift range 
lie on the fundamental plane determined by the rest of 
the quiescent population. This result suggests some com¬ 
monality of the origin of compact and more extended 
quiescent objects. Exploration of the links between en¬ 
vironment, abundance, and the Fundamental Plane may 
be further clues to the origin and evolution of compact 
galaxies. 

We thank the referee for careful reading of the 
manuscript and for suggesting the inclusion of Table 2. 
ID is supported by the Harvard College Observatory 
Menzel Fellowship and the Natural Sciences and Engi¬ 
neering Research Council of Canada Postdoctoral Fel¬ 
lowship (NSERC PDF-421224-2012). The Smithsonian 
Institution supports the research of MJG. HJZ gratefully 
acknowledges the generous support of the Clay Fellow¬ 
ship. 

REFERENCES 


Ahn, C. P., Alexandroff, R., Allende Prieto, C., et al. 2014, The 
Astrophysical Journal Supplement Series, 211, 17 
Arnouts, S., Cristiani, S., Moscardini, L., et al. 1999, Monthly 
Notices of the Royal Astronomical Society, 310, 540 
Barro, G., Faber, S. M., Perez-Gonzalez, P. G., et al. 2013, The 
Astrophysical Journal, 765, 104 
Barro, G., Trump, J. R., Koo, D. C., et al. 2014, The 
Astrophysical Journal, 795, 145 

Belli, S., Newman, A. B., & Ellis, R. S. 2014a, The Astrophysical 
Journal, 783, 117 

Belli, S., Newman, A. B., Ellis, R. S., & Konidaris, N. P. 2014b, 
The Astrophysical Journal, 788, L29 
Bezanson, R., van Dokkum, P. G., Tal, T., et al. 2009, The 
Astrophysical Journal, 697, 1290 
Bruzual, G., Chariot, S. 2003, Monthly Notices of the Royal 
Astronomical Society, 344, 1000 
Buitrago, F., Trujillo, I., Conselice, C. J., et al. 2008, The 
Astrophysical Journal, 687, L61 

Calvi, R., Poggianti, B. M., Vulcani, B. 2011, Monthly Notices 
of the Royal Astronomical Society, 416, 727 
Calzetti, D., Armus, L., Bohlin, R. C., et al. 2000, The 
Astrophysical Journal, 533, 682 
Carollo, C. M., Bschorr, T. J., Renzini, A., et al. 2013, The 
Astrophysical Journal, 773, 112 

Cassata, P., Giavalisco, M., & Guo, Y. 2011, The Astrophysical 
..., 743, 96 

Cassata, P., Giavalisco, M., Williams, C. C., et al. 2013, The 
Astrophysical Journal, 775, 106 


Cimatti, A., Cassata, P., Pozzetti, L., et al. 2008, Astronomy and 
Astrophysics, 482, 21 

Coil, A. L., Blanton, M. R., Buries, S. M., et al. 2011, The 
Astrophysical Journal, 741, 8 

Cool, R. J., Moustakas, J., Blanton, M. R., et al. 2013, The 
Astrophysical Journal, 767, 118 

Daddi, E., Renzini, a., Pirzkal, N., et al. 2005, The Astrophysical 
Journal, 626, 680 

Damjanov, I., Chilingarian, I., Hwang, H. S., Sz Geller, M. J. 

2013, The Astrophysical Journal, 775, L48 

Damjanov, I., Hwang, H. S., Geller, M. J., Sz Chilingarian, I. 

2014, The Astrophysical Journal, 793, 39 

Damjanov, I., McCarthy, P. J., Abraham, R. G., et al. 2009, The 
Astrophysical Journal, 695, 101 
Damjanov, L, Abraham, R., Glazebrook, K., et al. 2011, The 
Astrophysical Journal Letters, 739, L44 
Davari, R., Ho, L. C., Peng, C. Y., & Huang, S. 2014, The 
Astrophysical Journal, 787, 69 

Davies, L. J. M., Driver, S. P., Robotham, A. S. G., et al. 2014, 
Monthly Notices of the Royal Astronomical Society, 447, 1014 
Dawson, K. S., Schlegel, D. J., Ahn, C. P., et al. 2013, The 
Astronomical Journal, 145, 10 

Fasano, G., Marmo, C., Varela, J., et al. 2006, Astronomy and 
Astrophysics, 445, 805 

Ferre-Mateu, a., Vazdekis, a., Trujillo, I., et al. 2012, Monthly 
Notices of the Royal Astronomical Society, 423, 632 
Gehrels, N. 1986, The Astrophysical Journal, 303, 336 
Grogin, N. A., Kocevski, D. D., Faber, S. M., et al. 2011, The 
Astrophysical Journal Supplement Series, 197, 35 
Guzzo, L., Cassata, P., Finoguenov, A., et al. 2007, The 
Astrophysical Journal Supplement Series, 172, 254 
Ilbert, O., Arnouts, S., McCracken, H. J., et al. 2006, Astronomy 
and Astrophysics, 457, 841 

Ilbert, O., Salvato, M., Le Floc’h, E., et al. 2010, The 
Astrophysical Journal, 709, 644 
Ilbert, O., McCracken, H. J., Le Fevre, O., et al. 2013, 

Astronomy & Astrophysics, 556, A55 
Koekemoer, A. M., Aussel, H., Calzetti, D., et al. 2007, The 
Astrophysical Journal Supplement Series, 172, 196 
Koekemoer, A. M., Faber, S. M., Ferguson, H. C., et al. 2011, 

The Astrophysical Journal Supplement Series, 197, 36 
Kovac, K., Lilly, S. J., Cucciati, O., et al. 2010, The Astrophysical 
Journal, 708, 505 

Le Fevre, O., Vettolani, G., Garilli, B., et al. 2005, Astronomy 
and Astrophysics, 439, 845 

Lilly, S. J., Le Brun, V., Maier, C., et al. 2009, The Astrophysical 
Journal Supplement Series, 184, 218 
Longhetti, M., Saracco, P., Severgnini, P., et al. 2007, Monthly 
Notices of the Royal Astronomical Society, 374, 614 
Lopez-Sanjuan, C., Le Fevre, O., Ilbert, O., et al. 2012, 
Astronomy Sz Astrophysics, 548, A7 
Masters, K. L., Maraston, C., Nichol, R. C., et al. 2011, Monthly 
Notices of the Royal Astronomical Society, 418, 1055 
McCracken, H. J., Milvang-Jensen, B., Dunlop, J., et al. 2012, 
Astronomy & Astrophysics, 544, A156 
McLure, R. J., Pearce, H. J., Dunlop, J. S., et al. 2012, Monthly 
Notices of the Royal Astronomical Society, 428, 1088 
Muzzin, A., van Dokkum, P., Franx, M., et al. 2009, The 
Astrophysical Journal, 706, L188 
Nelson, E., van Dokkum, P., Franx, M., et al. 2014, Nature, 513, 
394 

Newman, A. B., Ellis, R. S., Bundy, K., Si Treu, T. 2012, The 
Astrophysical Journal, 746, 162 
Nipoti, C., Treu, T., Auger, M., Sz Bolton, A. 2009a, The 
Astrophysical Journal Letters, 706, L86 
Nipoti, C., Treu, T., Sz Bolton, A. S. 2009b, The Astrophysical 
Journal, 703, 1531 

Oogi, T., Sz Habe, A. 2013, Monthly Notices of the Royal 
Astronomical Society, 428, 641 

Peng, C. Y., Ho, L. C., Impey, C. D., Sz Rix, H.-W. 2010, The 
Astronomical Journal, 139, 2097 
Poggianti, B. M., Moretti, a., Calvi, R., et al. 2013a, The 
Astrophysical Journal, 777, 125 
Poggianti, B. M., Calvi, R., Bindoni, D., et al. 2013b, The 
Astrophysical Journal, 762, 77 
Raichoor, A., Mei, S., Stanford, S. A., et al. 2012, The 
Astrophysical Journal, 745, 130 


Quiescent Compact Galaxies at Intermediate Redshift in the COSMOS Field. The Number Density 


15 


Ross, N. P., Myers, A. D., Sheldon, E. S., et al. 2012, The 
Astrophysical Journal Supplement Series, 199, 3 
Santini, P., Ferguson, H. C., Fontana, a., et al. 2015, The 
Astrophysical Journal, 801, 97 

Saracco, P., Longhetti, M., & Gargiulo, A. 2010, Monthly Notices 
of the Royal Astronomical Society: Letters, 408, L21 
Sargent, M. T., Carollo, C. M., Lilly, S. J., et al. 2007, The 
Astrophysical Journal Supplement Series, 172, 434 
Scarlata, C., Carollo, C. M., Lilly, S., et al. 2007, The 
Astrophysical Journal Supplement Series, 172, 406 
Scoville, N., Aussel, H., Benson, A., et al. 2007, The 
Astrophysical Journal Supplement Series, 172, 150 
Scoville, N., Arnouts, S., Aussel, H., et al. 2013, The 
Astrophysical Journal Supplement Series, 206, 3 
Sersic, J. L. 1968, Atlas de galaxias australes (Cordoba, 

Argentina: Observatorio Astronomico) 

Shen, S., Mo, H., Sc White, S. 2003, Monthly Notices of the Royal 
Astronomical Society, 343, 978 

Simard, L., Trevor Mendel, J., Patton, D. R., Ellison, S. L., Sz 
McConnachie, A. W. 2011, The Astrophysical Journal 
Supplement Series, 196, 11 

Simard, L., Willmer, C. N. A., Vogt, N. P., et al. 2002, The 
Astrophysical Journal Supplement Series, 142, 1 
Stefanon, M., Marchesini, D., Rudnick, G. H., Brammer, G. B., Sz 
Whitaker, K. E. 2013, The Astrophysical Journal, 768, 92 
Szomoru, D., Franx, M., Sz van Dokkum, P. G. 2012, The 
Astrophysical Journal, 749, 121 


Taylor, E. N., Franx, M., Glazebrook, K., et al. 2010, The 
Astrophysical Journal, 720, 723 
Toft, S., van Dokkum, P., Franx, M., et al. 2007, The 
Astrophysical Journal, 671, 285 

Toft, S., Smolcic, V., Magnelli, B., et al. 2014, The Astrophysical 
Journal, 782, 68 

Trujillo, L, Cenarro, a. J., de Lorenzo-Caceres, A., et al. 2009, 
The Astrophysical Journal, 692, LI 18 
Trujillo, L, Conselice, C. J., Bundy, K., et al. 2007, Monthly 
Notices of the Royal Astronomical Society, 382, 109 
Valentinuzzi, T., Fritz, J., Poggianti, B. M., et al. 2010, The 
Astrophysical Journal, 712, 226 
van de Sande, J., Kriek, M., Franx, M., et al. 2013, The 
Astrophysical Journal, 771, 85 

van der Wei, a., Franx, M., van Dokkum, P. G., et al. 2014, The 
Astrophysical Journal, 788, 28 
van Dokkum, P. G., Franx, M., Kriek, M., et al. 2008, The 
Astrophysical Journal, 677, L5 

Wellons, S., Torrey, P., Ma, C.-P., et al. 2015, Monthly Notices of 
the Royal Astronomical Society, 449, 361 
Williams, C. C., Giavalisco, M., Cassata, P., et al. 2014, The 
Astrophysical Journal, 780, 1 

Zahid, H. J., Damjanov, L, Geller, M. J., Sz Chilingarian, 1. 2015, 
eprint arXiv:1501.04977 

Zirm, A. W., van der Wei, A., Franx, M., et al. 2007, The 
Astrophysical Journal, 656, 66 
Zolotov, A., Dekel, A., Mandelker, N., et al. 2014, eprint 
arXiv:1412.4783 



